Electronic structure of rare-earth impurities in GaAs and GaN 
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The electronic structures of substitutional rare-earth (RE) impurities in GaAs and cubic GaN 
are calculated. The total energy is evaluated with the self-interaction corrected local spin density 
approximation, by which several configurations of the open 4f shell of the rare-earth ion may be 
investigated. The defects are modelled by supercells of type REGa„_iAs n , for n = 4, 8 and 16. The 
preferred defect is the rare-earth substituting Ga, for which case the rare-earth valency in intrinsic 

material is found to be trivalent in all cases except Ce and Pr in GaN. The 3H > 2+ /-level is found 

above the theoretical conduction band edge in all cases and within the experimental gap only for 
Eu, Tm and Yb in GaAs and for Eu in GaN. The exchange interaction of the rare-earth impurity 
with the states at both the valence band maximum and the conduction band minimum is weak, 
one to two orders of magnitude smaller than that of Mn impurities. Hence the coupling strength is 
insufficient to allow for ferromagnetic ordering of dilute impurities, except at very low temperatures. 

PACS numbers: 75.50.Pp, 71.55.Eq, 71.70.Gm, 71.27. +a 
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I. INTRODUCTION 

The technology of electronics exploiting both charge- 
and spin-degrees of freedom of carriers is a rapidly de- 
veloping research field. 1,2 Diluted magnetic semiconduc- 
tors (DMS) are good candidates for new generations 
of fast, low dissipation, non- volatile integrated informa- 
tion and processing devices. In particular, ferromag- 
netic semiconductors offer a unique combination of mag- 
netic, semiconducting and optical properties, together 
with compatibility with existing microelectronics tech- 
nology. In the longer terms, the long electron life- 
times in semiconductors coupled with their strong inter- 
play with nuclear spins might be exploited for quantum 
computing* 3 " Several reports of ferromagnetically ordered 
DMS at room temperature (see e.g. Refs. and ref- 
erences therein) pose questions of fundamental impor- 
tance as to the nature of magnetic interactions in these 
materials, 6 ! 7 ! 8 ! 9 ! 10 . 11 : 12 

Most efforts have investigated Mn and other 3d im- 
purities in III-V and II- VI hosts, but recently also the 
discovery of high temperature ferromagnetism 13 ' 14 in di- 
luted Gd-doped GaN has gained considerable interest, 
implying that rare-earth (RE) ions may become a viable 
alternative to transition metals for spintronics applica- 
tions. In addition, rare-earth doping of wide-gap semi- 
conductors such as GaN and ZnO, is an active research 
area in itself due to the possibility of achieving tunable 
light emission^ arising from intra- /-shell optical transi- 
tions. 

The present work undertakes a systematic theoreti- 
cal investigation of the electronic structure of rare-earth 
dopants in GaAs and GaN, both in the zinc-blende crys- 
tal structure. Total energy calculations reveal Ga substi- 
tution as the most favorable rare-earth defect, with the 
rare-earth ion taking the iso-electronic trivalent config- 



uration in all cases studied except Ce and Pr in GaN. 
The effect of the exchange interaction of the rare-earth 
ion with the host states at the valence band maximum 
(VBM) and conduction band minimum (CBM) is calcu- 
lated, and is found surprisingly weak given the large spin 
moment formed on the rare-earth ion. A significant hy- 
bridization effect between the unoccupied /-states and 
the CBM states is observed. Furthermore, the position 
of the divalent acceptor level e(0/— ) is calculated and 
found outside the theoretical band gap, i.e. appearing 
as a resonance in the conduction bands. However, the 
e(0/— ) level falls inside the experimental band gap for 
Eu in GaN and for Eu, Tm and Yb in GaAs. 

The theoretical modelling of solids containing rare- 
earths poses severe problems due to the correlated na- 
ture of the / electrons in the incompletely filled 4/ shell. 
A description in terms of bands, as implemented in the 
conventional local spin density (LSD) approximationi&i 7 " 
invoked in the framework of density-functional theory^ 
leads to a large overestimation of the contribution of 
/-electrons to bonding and hence too small lattice con- 
stants. In addition, even if the band structures of LSD 
have no formal interpretation, they are often used as 
guidelines for what the elementary electron excitations of 
the solid are, but in the case of /-electron materials the 
approach is particularly poor. The / derived states ap- 
pear as a set of narrow bands situated close to the Fermi 
level, while experiments reveal their localized nature as 
atomic multiplet features over a broad energy ranged 

In the present work the electronic structure of rare- 
earth impurities in semiconductors is investigated with 
the self-interaction corrected (SIC) local spin density 
method, which constitutes an extension of the conven- 
tional LSD approximation by correcting for the spurious 
self-interaction of individual electrons^ For rare-earth 
atoms the self-interaction correction allows for a local- 
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ization of the /-electrons. In this way the atomic limit 
of the rare-earth / shell is implemented in conjunction 
with the itinerant limit of the other rare-earth valence 
electrons and the valence electrons of the semiconduct- 
ing host. In particular, the ground state configuration of 
the rare-earth /-shell is determined from comparison of 
the total energy corresponding to different valency sce- 
narios, as has already been shown in several applications 
to rare-earth solids, 21,22,23 including pressure induced va- 
lence transitionsi 24 i 25 i 26 

In section II of the present paper we outline the theo- 
retical aspects of the present work, notably the SIC-LSD 
total energy method as well as other calculational de- 
tails. In section III the results are presented for defect 
formation energies, magnetic interactions, and valency 
transition energies. Section IV contains the summary 
and conclusions of this work. 



II. THEORY 

A. The SIC-LSD total energy method 

The total energy functional of the LSD approxima- 
tion is renowned for its chemical accuracy in describing 
conventional solids with weakly correlated elect rons>iSiii 
The self-interaction correction is included to facilitate an 
accurate description of the localized / electrons of rare- 
earths. Specifically, the SIC-LSD 20 total energy func- 
tional is obtained from the LSD as: 



^SIC-LSD 

AEgic 
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occ. 



a 

occ. 

Ec 



(i) 

(2) 
(3) 



where a labels the occupied states and S^ IC is the self- 
interaction correction for state a. As usual, E LSD can 
be decomposed into a kinetic energy, T, a Hartree en- 
ergy, U, the interaction energy with the atomic ions, 
V ex t, and the exchange and correlation energy, E xc *£ 
The self-interaction is defined as the sum of the Hartree 
interaction and the exchange-correlation energy for the 
charge density of state a: 



rSIC 



U[n a ] + E xc [n a }. 



(4) 



For itinerant states, <5f /c vanishes identically, while for 
localized (atomic-like) states the self-interaction may be 
appreciable. This correction constitutes a negative en- 
ergy contribution for an /-electron to localize, which 
then competes with the band formation energy gained 
by the /-electron if allowed to delocalize and hybridize 
with the available conduction states. In rare-earths the 
self-interaction correction ranges from S a ~ 0.8 eV per 



/-electron in Ce to S a ~ 1.5 eV per /-electron in Yb, re- 
flecting the contraction of the /-orbitals through the se- 
ries. The volume dependence of 5 a is rather weak, hence 
the overbinding of the LSD approximation for narrow 
band states is reduced when localization is allowed. The 
last term in Eq. (1) is the spin-orbit energy, where for 
each occupied state a: 



(•5) 



is added to the energy functional. We employ the atomic 
spheres approximation, whereby the crystal volume is di- 
vided into slightly overlapping atom-centered spheres of 
a total volume equal to the actual volume. In (JjjJ, the 
angular momentum operator, I — r x p, is defined in- 
side each atomic sphere, with r given as the position 
vector from the sphere center. The spin-orbit interac- 
tion couples the band Hamiltonian for the spin-up and 
spin-down channels, i.e. a doubled secular problem must 
be solved. Other relativistic effects are automatically in- 
cluded by solving the scalar-relativistic radial equation 
inside spheres. The spin-orbit parameter, 

£( r ) = ^1^- 

c 2 dr 

in atomic Rydberg units, is calculated from the self- 
consistent potential. The SIC-LSD energy functional in 
Eq. (1) appears to be a functional of all the one-electron 
orbitals, but can in fact be viewed as a functional of the 
total (spin) density alone, as discussed in Ref. lU 

The advantage of the SIC-LSD energy functional is 
that it allows for different valency scenarios to be ex- 
plored. By assuming atomic configurations with differ- 
ent total numbers of localized states, self-consistent min- 
imization of the total energy leads to different local min- 
ima of the same functional, E SIC ~ LSD in Eq. (1), and 
hence their total energies may be compared. The config- 
uration with the lowest energy defines the ground state 
configuration. Note, that if no localized states are as- 
sumed, E SIC ~ LSD coincides with the conventional LSD 
functional, i.e., the Kohn-Sham minimum of the E LSD 
functional is also a local minimum of E SIC ~ LSD . The 
interesting question is, whether competing minima with 
a finite number of localized states exist. This is usually 
the case in /-electron systems^ and some 3d transition 
metal compounds^ where the respective / and d orbitals 
are sufficiently confined in space to benefit appreciably 
from the SIC. 

The SIC-LSD still considers the electronic structure of 
the solid to be built from individual one-electron states, 
but offers an alternative description to the Bloch picture, 
namely in terms of periodic arrays of localized atom- 
centered states (i.e., the Heitler-London picture in terms 
of Wannier orbitals). Nevertheless, there still exist states 
which will never benefit from the SIC. These states retain 
their itinerant character of the Bloch form, and move in 
the effective LSD potential. This is the case for the host 
valence band states in the systems studied here. The 
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SlC-term only affects the rare-earth ions. The resulting 
many-electron wavefunction will consist of both localized 
and itinerant states. The nominal valency of a rare-earth 
ion is defined as the number of valence electrons available 
for band formation: 

N va i — Z — N core - NsiC, 

where Z is the nuclear charge, N core the number of core 
electrons, and N$ic the number of self-interaction cor- 
rected /-electrons. 

In contrast to the LSD Kohn-Sham equations, the SIC 
electron states, which minimize E , experience 

different effective potentials. This implies that to mini- 
mize E SIC ~ LSD , it is necessary to explicitly ensure the 
orthonormality of the one-electron wavefunctions by in- 
troducing a Lagrangian multipliers matrix. Furthermore, 
the total energy is not anymore invariant with respect to 
a unitary transformation of the one-electron wavefunc- 
tions. Both of these aspects make the energy minimiza- 
tion more demanding to accomplish than in the LSD 
case. The electron wavefunctions are expanded in linear- 
muffin-tin-orbital (LMTO) basis functions, 2 - and the en- 
ergy minimization problem becomes a non-linear opti- 
mization problem in the expansion coefficients. Further 
details of the present implementation can be found in 
Ref. S 



B. Calculational Details 

Alloys of GaAs (or GaN) and rare-earth ions of 
the type RE^Gai-^As (RE^Gai-xN) and RE^GaAsi-j; 
(RE^GaN^^) were considered by supercell modelling, 
and total energies evaluated for x = 1/4, 1/8 and 1/16. 
The experimental lattice constants of pure GaAs and 
GaN in the zinc-blende structure were used (a — 5.65 
A for GaAs, a = 4.45 A for GaN), i.e. the effects of 
lattice expansion with doping were assumed negligible. 
The atomic positions were kept fixed in all calculations, 
except for a few test cases discussed in the next section, 
i.e. no systematic study of the effects of relaxations of 
nearest neighbours of the impurities were attempted. Ex- 
perimentally, for Er-doped GaN the Er-N distance has 
been found ~ 10% larger than the pure host Ga-N bond 
length piii 

The itinerant states were sampled using 95, 22 and 14 
k-points in the irreducible wedge of the Brillouin zone 
of the x — 1/4, 1/8 and 1/16 supercell, respectively. In 
GaAs, the id semicore states of Ga and As were treated 
selfconsistently in a separate energy panel, which also in- 
cluded the 5s and 5p semicore states of the rare-earth. In 
the case of the GaN host, the best description is obtained 
with only the 5s semicore states of the rare-earth impu- 
rity treated in the separate semicore panel. Due to the 
smaller lattice constant, there is a pronounced interaction 
between the N 2s, 2p states and the Ga 3d and rare-earth 
5p states, for which reason the latter were treated as va- 



lence electrons to describe properly their hybridization 
into the GaN valence bands. 



III. RESULTS AND DISCUSSION 
A. Cohesive Properties 

Figures n an d 13 show the energy differences calculated 
for trivalent, divalent and tetravalent configurations of 
the rare-earth impurities, in GaAs and GaN, respectively, 
at concentration x = 1/16. The rare-earth ion is taken 
to substitute for Ga. The preferred valency is seen to be 
trivalent for all rare-earth impurities in GaAs, and for all 
but Ce and Pr in GaN, in good accord with the triva- 
lency of the Ga atom which the rare-earth replaces. In 
particular, the Eu and Yb impurities remain trivalent (by 
a margin of more than 0.5 eV) in spite of the tendency of 
these elements to become divalent, e.g. in the elemental 
metallic phase. Luminescence experiments show that Eu 
in GaN is trivalent^ although there have been indica- 
tions of divalent Eu in the surface region of Eu-doped 
GaN^ 

In GaAs, the valence energy difference is almost inde- 
pendent of concentration, except for Ce, for which case, 
at concentrations x = 1/4 and x = 1/8 the preferred con- 
figuration is f° (tetravalent). This means that /-band 
formation is favored for Ce impurities at this large con- 
centration, i.e. instead of having one localized /-electron 
on each Ce atom, a narrow /-band is formed close to 
the conduction band minimum, and filled with one elec- 
tron leading to a metallic alloy. However, in the most 
diluted case studied here, x = 1/16, the localization of 
one /-electron on each Ce is more favorable, leading to 
a ground state with cerium ions in the trivalent f 1 con- 
figuration and a semiconducting alloy, which we hence 
assume to prevail also at Ce concentrations of even lower 
value. 

In GaN both Ce and Pr are found to be tetravalent, 
even at x — 1/16, i.e. again implying the formation of 
narrow /-bands straddling the conduction band edge and 
metallicity. The band formation is favored by the small 
host lattice constant, which particularly affects the least 
tightly bound /-electrons in the early lanthanidcs. Pre- 
sumably, further dilution beyond what has been studied 
here, will eventually lead to the formation of localized 
trivalent ions also for Ce and Pr in GaN. However, also 
effects of relaxations of the nearby host atoms are impor- 
tant. As a check, a similar valency energy difference was 
calculated for Ce in GaN with the nearest neighbor N 
atoms moved 10% of the bond length away from Ce, and 
indeed the trivalent configuration in this case had a lower 
energy (by 0.14 eV) than the tetravalent configuration. 
The detailed analysis of the localization-delocalization 
transition upon dilution would need to incorporate be- 
sides the relaxation of atomic positions in the vicinity 
of the impurity, also issues of impurity disorder and/or 
clustering, which is not attempted here. 
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Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 




Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 



FIG. 1: (Color online) The energy difference between tetrava- 
lent and trivalent rare-earth impurities substituting for Ga in 
GaAs (circles), and for selected cases between divalent and 
trivalent ions (triangles). The calculations were performed 
for supercells corresponding to rare-earth dilution x = 1/16. 
The positive sign implies that the trivalent configuration is 
the ground state. 



FIG. 2: (Color online) The energy difference between tetrava- 
lent and trivalent rare-earth impurities substituting for Ga 
in cubic GaN (circles), and for selected cases between di- 
valent and trivalent ions (triangles). The calculations were 
performed for supercells corresponding to rare-earth dilution 
x = 1/16. 



The trends of the formation energy of charge neutral 
substitutional rare-earth ions in GaAs and GaN are de- 
picted in Fig. |21 The formation energy is defined relative 
to the pure solids GaAs (GaN) and metallic rare-earth: 

Eform = -E(REGai5Asi6) + flCa — 

-E(Gai 6 Asi 6 ) - hre, (6) 

where the two energies, E^REGaisAsig), and 
£J(Gai6Asi6), are total energies calculated for the 
REGai5Asi6 and GaigAsi6 supercells. The chemical 
potential for the rare-earth is taken as the total energy 
for the elemental metallic phase /j,re = E(RE), likewise 
calculated with the SIC-LSD methodiSi /iQa is the 
chemical potential for Ga, which depends on the growth 
conditions of the doped semiconductor. During growth 
the chemical potentials for Ga and As are constrained by 
fJ-Ga + ^As = .E(GaAs). For Ga-rich conditions, excess 
Ga is assumed to be in the metallic phase (orthorhombic 
a-Ga, Cmca symmetry), i.e. fica = S(a-Ga), while for 
As-rich (N-rich) conditions excess As (N) is assumed 
to exist as solid As (rhombohedral a- As, symmetry 
R3m) and N2 gas, respectively, i.e. has = E(a-As) 
and yitjv = ^E{N2). The formation energies in Fig. [31 
refer to the pnictide-rich growth conditions. For Ga-rich 
conditions, the formation energies of Fig. should be 
raised by 1.2 eV in the GaAs host and by 2.0 eV in 
the GaN host, respectively. The reference state of the 
rare-earth metal is trivalent, except for Eu and Yb, 
which are divalent, due to the advantage of attaining 
the half-filled f 7 configuration, and full f 14 shell, 
respectively. 21 This causes a relatively high formation 



energy of Eu(/ 6 ) defects in the GaAs and GaN hosts. 
Apart from the anomaly of Eu, the formation energies 
are seen to exhibit a monotonously decreasing trend 
through the rare-earth series with a shallow minimim 
towards the later rare-earths Dy-Tm. The relatively low 
formation energies are consistent with the high level of 
doping (~ 1%), which may be obtained for rare-earths 
in GaN. 15 The formation energies will be somewhat 
lowered by relaxations of the host atomic positions in the 
vicinity of the impurity, which were not considered here. 
A check of the relaxation effects was performed for Gd 
in GaN using the LSD approximation as implemented in 
the full-potential version of the LMTO approach of Ref. 
liiil This calculation revealed an 1 1 % outward relaxation 
of the nearest neighbor N atoms of the Gd impurity, 
which is in excellent agreement with the experimental 
investigations of Er ions in GaN, finding 10 % outward 
relaxation^ The use of the LSD approximation for this 
particular system is justified by the similarity of the 
electronic structure of Gd in LSD and SIC-LSD (due 
to the half- filled /-shell), but cannot be relied upon in 
general. 

A few cases of rare-earth ions substituting for As in 
GaAs have been considered: Pr, Eu and Tb. For these, 
the formation energies are substantially higher, namely 
9.0 eV, 8.4 eV, and 8.0 eV, respectively, again for As-rich 
conditions (the values are 1.2 eV lower for Ga-rich con- 
ditions). We conclude, that indeed the As-substitution 
is unfavorable compared to Ga-substitution, as also evi- 
denced by experiments. 15 

In Fig. 0| the electronic density of states is illustrated 
for the case of Eu in GaN. The calculated gap is 2.0 eV 
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Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 

FIG. 3: (Color online) Formation energies (cf. Eq. ©) for 
substitutional rare-earth ions in GaAs (circles) and GaN (as- 
terix) under pnictide-rich conditions (see text for discussions) . 
For Ga-rich conditions, the formation energies are larger by 
1.2 eV (GaAs) and 2.0 eV (GaN). 



(in the presence of Eu) , defined as the separation of the 
host Tg valence band maximum state and the Tg con- 
duction band minimum state. The figure reveals a sharp 
resonance inside the gap region originating from the sin- 
gle Eu majority /-state, which is not localized in the / 6 
configuration, and hence appears as an unoccupied impu- 
rity band state. We discuss this state later, in subsection 
III.C. The minority Eu /-states are seen to hybridize into 
the conduction bands 4-5 eV above the VBM. The en- 
ergy position of the localized states are not shown (see 
discussion in subsection fff.D). 



B. Magnetic Properties 

The magnetic interaction of the rare-earth dopants 
with the host band states is of crucial importance for the 
potential spintronics applications. We analyze this by 
monitoring the exchange splitting, AE = e(j) — e(T), of 
the band states at the valence band maximum (VBM) 
and the conduction band minimum (CBM), both of 
which occur at the T point in the pure hosts as well as for 
the x = 1/16 impurity supercells. Here the spin-direction 
| is taken as the direction of the rare-earth spin (majority 
spin direction). At the VBM, the topmost j = 3/2 four- 
fold degenerate valence state Tg splits into four distinct 
energies separated by a few meV. Two of these are es- 
sentially pure spin-up and spin-down, respectively, (cor- 
responding to the rrij = ±3/2 quantum numbers), and 
Fig. [S] shows the calculated differences in energy between 
these two states. The CBM state is the twofold r 6 state, 
which is antibonding s-like, and its splitting in the pres- 




E (eV) 

FIG. 4: Electron density of states (DOS) for Eu in GaN 
(x = 1/16). Full line shows the total DOS of the x = 1/16 
supercell in units of states per eV and supercell. The dashed 
line shows the DOS projected onto the Eu site. The zero of 
energy is placed at the VBM. 




Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 



FIG. 5: (Color online) Spin-splitting (in meV) of the electron 
bands at the valence band maximum for rare-earth dopants in 
GaN (circles) and GaAs (diamonds) for x = 1/16. A positive 
sign implies that the band state with its spin aligned with the 
rare-earth spin has a lower energy than the state with its spin 
anti-parallel to that of the rare-earth ion. The sign changes 
of the spin-splitting for GaAs at Ho is caused by the mixing 
in of the unoccupied rare-earth impurity /-bands into the top 
of the valence band. 



ence of magnetic rare-earth ions is depicted in Fig. 

The spin-splitting of the top of the valence band in 
GaN is seen to peak around Nd at a value of 45 meV 
(Fig. |5J), and falling smoothly towards a value of 6 meV 
in Yb. In contrast, the spin-splitting of the valence band 
of GaAs fluctuates more, which is due to the interaction 



6 



150 




-100- 



Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 

FIG. 6: (Color online) Spin-splitting (in meV) of the elec- 
tron bands at the conduction band minimum for rare-earth 
dopants in GaN (circles) and GaAs (diamonds) for x = 1/16. 
A positive sign implies that the band state with its spin 
aligned with the rare-earth spin has a lower energy than the 
state with its spin anti-parallel to that of the rare-earth ion. 
The sign changes of the spin-splitting at the conduction band 
minimum around Sm and Ho is caused by the unoccupied 
rare-earth impurity /-bands crossing the conduction edge (cf. 
Fig. 0. Note the larger energy scale compared to Fig. |5] 



with the unoccupied /-states. The latter appear in the 
SIC-LSD description as impurity bands, moving down 
with respect to the host bands, as one moves across the 
first half of the lanthanide series, and again similarly for 
the second half of the series. This is illustrated for GaN 
in Fig. [7| The hybridisation effect of the /-bands with 
the VBM p-bands, as well as with the CBM s-bands can 
be substantial^ if the /-bands are close in energy to 
the gap edges. In the case of GaAs the /-impurity band 
actually merges with the VBM for Sm and Eu dopants, 
while for Pm dopants the /-band is just above the VBM. 
Due to level repulsion the majority spin p-band at the 
VBM for GaAs:Pm is pushed significantly more down. 
Similarly, for the other spin channel in the second half 
of the lanthanide series, the /-band crosses the VBM 
between Ho and Er, and for Yb the /-like eigenvalue 
at the T-point is below the 4-fold degenerate host VBM 
states (but only at T - the impurity band is about one 
third filled at self-consistency). 

The sign of the spin-splitting at the top of the va- 
lence band is in all cases, except Ho in GaAs, positive, 
i.e. the VBM p-states tend to align with the rare-earth 
spin. Holes will therefore predominantly be created in 
the minority spin channel, i.e. the p-type carriers will be 
spin- aligned with the rare-earth dopant. A similar small 
interaction was found for Gd in GaN in Ref. |3f| which 
did not include spin-orbit interaction. However, in that 
work the opposite sign of the spin-splitting was reported. 
The origin of this sign difference most probably lies in 



the different treatment of the occupied Gd /-states. In 
Ref. |35|the /-states are treated as band states, which fall 
inside the occupied valence bands and interact strongly 
with the valence band top states. The repulsion effect on 
the majority spin states at the VBM is so strong, that 
these are pushed above the minority VBM state, leading 
to a sign change of AE. 

In contrast to the rather modest spin-splitting ob- 
tained for the rare-earth dopants in the present work, 
for Mn in GaN and GaAs it is much stronger. The 
spin-splitting induced by the Mn localized spins has been 
calculated^ by the SIC-LSD method to be AE = -0.55 
eV and -0.43 eV in GaN and GaAs (for x = 1/16), re- 
spectively, i.e. one or two orders of magnitude larger 
(and of opposite sign) than the splittings shown in Fig. 
13 For Mn in GaAs and GaN, the overlap of the Mn d- 
states (in particular those of ti symmetry) with the host 
states at the VBM is significantly larger than that of the 
rare-earth ions. Hence, for a Mn impurity with localized 
bonding and Mn d dominated states, the highest state at 
the VBM becomes the antibonding host-^2-Mn d band- 
like resonance, which is pushed up. At the same time, in 
the minority channel a bonding host-£2-Mn d resonance is 
formed and moving down in energy, altogether leading to 
a negative spin-splittin g A E = e(l) — e(T) 1 ^ The calcu- 
lations reported in Ref. [ll] did not include the spin-orbit 
interaction, but the conclusions of that work prevail when 
spin-orbit interaction is included^ In the case of rare- 
earth dopants the / orbitals are much better confined 
to the rare-earth site and their hybridizations are signif- 
icantly weaker, so that the spin-splitting becomes dom- 
inated by the direct exchange between host-states and 
rare-earth site (described through the spin-dependent ef- 
fective potential). The only exception occurs when the 
unoccupied /-bands accidentally fall close to the host 
band edges, as discussed. Hence we conclude that for a 
p-type GaAs or GaN host there will be a much weaker 
spin polarization effect induced on the hole carriers by 
rare-earth dopants than when transition metal dopants 
are used. 

GaN is most often n-type, so it is more appropriate to 
investigate the rare-earth induced effects on the conduc- 
tion band states. This is a more difficult matter, since 
the LDA bands of the unoccupied conduction states are 
known to be inaccurate. Most importantly, the funda- 
mental gaps are seriously underestimated, but in some 
cases the entire shape of the conduction bands is in er- 
ror, e.g. in GaAsi^l In the present approximation, the 
GaAs and GaN band gaps are calculated to be 0.11 eV 
and 1.77 eV, respectively, as compared to the experi- 
mental values of 1.42 eV— and around 3.3 eV—, and to 
full-potential LMTO values of 0.21 eVi& (see also Ref. 
l4l|) and 2.2 eV~ A further complication is the appear- 
ance of the unoccupied /-bands of the rare-earth dopant, 
which interact with the host conduction band states, in 
particular with the s-like CBM. 35 The physical meaning 
of these /-bands depends on whether the band picture 
applies for an /-like electron added onto the rare-earth 
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Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 

FIG. 7: (Color online) Trends in energy band structure 
around the fundamental gap for RE^Gai-^N, x — 1/16, su- 
percells. The VBM is at zero energy, full circles mark the 
unoccupied /-bands, while the squares and dotted line mark 
the calculated CBM. The /-band is defined by the smallest 
and largest /-dominated eigenvalue at the F-point (consider- 
ing majority spin only, for first half of the rare-earth-series), 
and its width is mainly due to spin-orbit splitting. The ex- 
perimental gap (3.3 eV) is marked with a dashed line. The 
calculated gap (for pure GaN) is 1.77 eV. The diamonds mark 
the calculated positions of the e(0/— ) levels in the localized 
picture, cf. Table I and subsection IIFC. 



ion, which is doubtful. For example, the observation of 
distinct multiplct lines in the photoluminesccncc spectra 
of rare-earth doped GaNi£ clearly demonstrates the lo- 
calized nature of the /-electrons. See also the discussion 
in the next subsection. 

The spin-splitting of the CBM states seen in Fig. El 
is minute in GaAs, but much stronger in GaN, in par- 
ticular where the unoccupied /-bands cross the CBM, 
which occurs around Sm and Ho (see Fig. [7J. In the 
mid-lanthanide series, for Eu, Gd, Tb, Dy, and Ho, the 
exchange interaction is antiferromagnetic, i.e. the minor- 
ity spin CBM state is lower than the majority spin state. 
This is again caused by the level repulsion effect. In the 
case of Eu, the unoccupied band of the seventh majority 
spin electron has moved below the conduction band edge 
thus pushing the majority spin s-state at CBM up. For 
Gd, Tb, Dy and Ho, the unoccupied /-band states are of 
minority spin, and situated above the CBM, thus push- 
ing the minority spin s-like CBM state down below the 
majority spin s-like CBM state. This behavior was also 
demonstrated for Gd in hypothetical zinc-blende GdN in 
Ref. |H 

To further analyze the exchange interaction of the rare- 
earth ion with the p-like VBM, the relevant exchange 



Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 

FIG. 8: (Color online) Exchange interaction parameter N/3 
of the valence band maximum for rare-earth dopants in GaN 
(circles) and GaAs (diamonds) for x = 1/16. 



integral, J p f is extracted using the expression^ 

AEvbm 



N/3 = NJ pf 



2Sx 



(7) 



Similarly, the exchange-splitting of the CBM s-like state 
can be quantified through an s — / interaction parameter: 



Na = NJ sf 



AE, 



CBM 



2Sx 



(8) 



In these expressions, S is the total spin of the rare-earth 
ion, x the dopant concentration, and N the number of 
host formula units within the normalization volume en- 
tering the definition of exchange integrals^ Figures E| 
and El show the values of these parameters through the 
rare-earth series. 

The exchange integral parameters are seen to reflect 
the behaviors already discussed for the spin-splittings 
in Figs. El and El For Gd in GaN we may compare 
to the exchange parameters extracted in Ref. |35| for 
GdN, whose estimated coupling parameters N(3 — —0.01 
eV and Na = —0.4 eV are showing somewhat weaker 
coupling to the VBM but significantly stronger coupling 
to the CBM, as compared to our respective values of 
N(3 — 53 meV and JVq = —46 meV, calculated for the 
diluted Gd in GaN (see Figs. Eland 03). 

In contrast to the rather modest interaction parame- 
ters obtained for the rare-earth dopants here, for dilute 
Mn in GaN and GaAs, the exchange interaction of the Mn 
localized spins with the VBM was obtained as N(3 = —2.5 
eV in GaN and N0 = -1.5 eV in GaAs4i 

In conclusion, it appears very unlikely that isolated 
rare-earth dopants by themselves can cause the room 
temperature ferromagnetism, let alone the giant magne- 
tization amplifications, observed in the experiments by 
Ref. E3 
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Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 



FIG. 9: (Color online) Exchange interaction parameter Na 
of the conduction band minimum for rare-earth dopants in 
GaN (circles) and GaAs (diamonds) for x = 1/16. 



C. Acceptor levels 

The behavior elucidated in Figs. and OH can only 
be considered a trend study, since the absolute posi- 
tions of the conduction bands are not accurate. The 
LSD is known to substantially underestimate the band 
gap of semiconductors^ In the SIC-LSD approach the 
host band states are described by the LSD approxima- 
tion, however applied to the self-consistent charge den- 
sity as obtained by the minimization of the SIC-LSD en- 
ergy functional, Eq. ^ which deviates somewhat from 
the self-consistent charge density obtained by minimizing 
the LSD total energy functional. Generally, the SIC-LSD 
leads to larger magnetic moments and larger /-electron 
occupation than the LSD, although sometimes the differ- 
ence can be rather small. Hence we may anticipate that 
the gap problem of LSD prevails in SIC-LSD. 

Even more problematic is the description of the unoc- 
cupied /-states as weakly correlated LSD bands, which 
is in conflict with the localized picture used for the oc- 
cupied part of the /-manifold. The LSD potential is de- 
rived from the limit of nearly-free electrons, i.e. fast elec- 
trons moving in an effective potential as given by the LSD 
approximation. This is a good approximation for most 
cases, e.g. for the occupied states of the host materials, 
which are built from extended s- and p-orbitals of the 
constituent atoms. The approximation is less adequate 
for the spatially confined /-orbitals of the rare-earth 
dopant, for which the overlap with neighbouring impu- 
rity atoms is small. Hence, an extra added /-electron 
would tend to stay for a long time on one rare-earth ion 
before hopping to a neighbour. The local environment 
will have time to adjust to the presence of the extra elec- 
tron, and the effective potential will be modified. The 
same effect is well known for localized gap states in semi- 



conductors, for which the transition state approximation 
is often invoked, i.e. the defect energy level is calculated 
by occupying the defect state with half of an electron. 
In this subsection we will investigate various approaches 
to the calculation of the /-related acceptor levels e(0/— ) 
corresponding to the addition of one /-electron to the 
/-shell. 

In the delocalized limit, as discussed, the added /- 
electron is put in the lowest unoccupied /-impurity band, 
which is positioned at 

(-LSD —< fn+l\HLSD\fn+l >, (9) 

where \f n +i > is the one-particle /-state, and Hlsd is 
the LSD Hamiltonian evaluated with the self-consistent 
charge density as obtained by minimization of the 
fiSic-LSD f unc tional. Generally, several unoccupied /- 
states exist, and the lowest value of the above matrix 
element must be taken to represent the minimal addition 
energy. 

In the opposite limit of localized / states one may cal- 
culate the position of the e(0/— ) level of a divalent rare- 
earth ion with one additional localized /-electron, by an 
SIC-LSD total energy difference: 

e(0/-) = E(r +1 ,-)-E(f n ,0) (10) 

Here, E(f n ,0) is the total energy of the charge neutral 
trivalent rare-earth ion, while E(f n+1 , — ) refers to the to- 
tal energy of a divalent rare-earth ion with an extra elec- 
tron in the supercell, and a homogeneous positive back- 
ground charge of +e added. In the approximation im- 
plied by Eq. i|10|) it is thus assumed that the presence of 
an extra electron on the rare-earth ion may be artificially 
compensated by the positive background charge, render- 
ing the supercell altogether electrically neutral, which is 
necessary in order to have a well defined total energy 
minimization problem. Obviously, spurious electrostatic 
interactions are introduced by this approach, and the re- 
liability of the scheme must be checked. The same ap- 
proach is often invoked for other charged impurities in 
semiconductors 

A good test of the approximation in Eq. I|1(J|I is pro- 
vided by the charge neutral SIC-LSD approximation for 
the / addition level. The latter is defined as the SIC- 
LSD total energy difference between divalent and triva- 
lent rare-earth ions in a charge neutral supercell, i.e. the 
extra electron added to the rare-earth is simply supplied 
from the highest valence band state. This corresponds to 
the energy differences displayed by the E(II) — E(III) 
curve in Figs. ^and^ 

e(0/-) apP r - evBM = E(f n+1 ,0) - E(f n , 0). (11) 

In this approximation the net accumulation of negative 
charge in the vicinity of the rare-earth ion is simply 
supplied from the Fermi level, however due to the fi- 
nite supercell, uncontrollable interactions between the 
negative /-ion and the hole state introduced will occur. 
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The charge density of the host hole state does not corre- 
spond to a completely homogeneous charge density, and 
furthermore, the missing electron charge influences also 
the exchange-correlation part of the energy functional. 
Hence the two approaches differ, but as we will show give 
quite similar results, so that we can conclude that indeed 
the Coulombic interaction caused by the finite supercell 
is small. 

Finally, we can also improve on the delocalized descrip- 
tion of the / acceptor level by considering the transition 
state. The position of the lowest unoccupied /-band at 
half occupancy may be approximated as the average of 
the position of this band at zero and full occupancy: 

4s = l^LSD(II) + e LS D(III)). (12) 

At zero occupancy, €lsd{HI) corresponds to Eq. J§J, 
and at full occupancy €lsd(H) corresponds to the same 
expression, however evaluated with the LSD Hamiltonian 
for the self-consistent charge density corresponding to 
the divalent rare-earth configuration. In this expression, 
|/n+i > must be that particular /-state, which is delo- 
calized in the trivalent but localized in the divalent con- 
figuration. The difference between the two LSD eigenval- 
ues entering the average can be identified as an effective 
Coulomb U parameter, and the transition state energy is 
therefore equivalent to a shift by of the LSD unoccu- 
pied band. 



Host / gap 


dopant 


<o/-) 


^(0/ )appr 


tLSD 






Pm 


2.30 


2.13 


0.37 


2.39 


GaAs 


Sm 


1.47 


1.28 


0.13 


1.96 




Eu 


0.89 


0.80 


0.00 


1.88 


expt: 


Gd 


4.17 


4.20 


1.71 


4.98 


1.42 


Tb 


3.25 


3.23 


0.61 


3.13 




Ho 


2.13 


2.15 


0.24 


2.83 


theo: 


Er 


1.85 


1.70 


0.00 


2.59 


0.11 


Tm 


1.12 


1.20 


0.00 


2.41 




Yb 


0.89 


0.97 


0.00 


2.38 




Nd 


4.21 


4.43 


2.91 


5.09 


GaN 


Pm 


4.03 


4.06 


2.56 


4.90 




Sm 


3.38 


3.11 


2.12 


3.69 


expt: 


Eu 


2.44 


2.15 


0.61 


2.79 


3.3 


Ho 


4.62 


4.44 


1.99 


4.57 




Er 


4.32 


4.19 


1.65 


4.44 


theo: 


Tm 


3.99 


3.89 


1.17 


4.37 


1.77 


Yb 


3.62 


3.45 


0.57 


3.78 



TABLE I: Calculated defect levels of divalent rare-earth 
dopants in GaAs and GaN. First column gives the host with 
experimental and theoretical gap values, second column the 
dopant, third column the calculated defect level on the basis 
of the SIC-LSD total energy difference, Eq. IjlOjl . fourth col- 
umn the same with the approximate expression, Eq. Hill , i.e. 
neglecting the charging of the near defect region, fifth column 
the LSD position of the lowest unoccupied /-band, Eq. @, 
and sixth coulmn the transition state as defined in Eq. 1121 . 
All numbers are in eV and relative to the VBM. 

In Table I we compare for selected cases the /-related 



e(0/— ) level in Eq. I|10|) with the approximate expres- 
sions given by equations © , !fl"T|) and (fF2*)l . The theoret- 
ical and experimental host gaps are quoted for compar- 
ison. The most consistent interpretation of the present 
results is to compare the calculated levels with the the- 
oretical band gap, but often comparisons to the ex- 
perimental gap in fact are done. After all, the band 
gap discontinuity** problem persists in the expressions 
CHll-G2J- One observes that in all cases the negatively 
charged rare-earth impurity has a higher energy than the 
theoretical band gap, and hence appears as a scattering 
resonance in the conduction bands. The implication is 
that even in ro-type material the rare-earth ions will re- 
main trivalent. If compared to the experimental band 
gap, only for Eu, Tm and Yb in GaAs and for Eu in GaN 
does the e(0/— ) level fall below the experimental con- 
duction band edge, while for Sm the e(0/— ) level comes 
rather close to the experimental conduction edge. Photo- 
luminescence of (wurtzite) GaN:Eu thin films reveals no 
signature of divalent Eu ions^ thus corroborating the 
conclusions drawn by comparison to the theoretical gap. 

Through the rare-earth series, the approximate expres- 
sion in Eq. {TT} differs from Eq. {TUJl by only 0.1 - 0.3 
eV, which is reassuring, since it shows that the spurious 
charge interaction effects of the divalent rare-earth ion in 
the finite supercell are quite small. The correction with 
respect to the LSD band position is significant, generally 
several eV, i.e. the local Coulombic corrections due to 
the presence of the added localized /-electron are sub- 
stantial. As explained, the transition state construction 
seeks to take this into account, and indeed brings the 
level position in much better agreement with the expres- 
sion in Eq. (|10|) . although it tends to mostly overshoot, 
in some cases by more than 1 eV. 



D. f removal energies 

The calculation of the energy position of the localized 
states in SIC-LSD theory has always been a matter of 
concern. The /-removal energies one observes in pho- 
toemission include all the atomic multiplets of the / n_1 
ion left behind, which can lead to a spectral function 
with f- related features over a ~ 10 eV range, see e.g. 
Ref. |43. Even calculating the energy of the most stable 
state with one /-electron removed is difficult. First of 
all, this will usually not be the ground state for an N — 1 
electron system, so applying a density functional theory 
based formalism poses a problem. But also in practice, 
the self-consistency process for a rare-earth with a hole 
created in the / shell leads to unwanted effects, since 
the enhanced Coulomb attraction due to the missing /- 
electron will pull down the unoccupied /-levels below the 
Fermi level. Hence, if the initial state corresponds to a 
localized /" configuration, the effective rare-earth config- 
uration in the final state will not be / n_1 , but somthing 
in between /™ _1 and /", often closer to the latter, and 
hence a too small removal energy is obtained (photoemis- 
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sion experiments do interprete certain features in terms 
of such 'well screened' peaks^). On the other hand, just 
knocking out the /-electron but not driving the charge 
density to self-consistency overestimates the / removal 
energy by not allowing for screening of the charge hole 
by non-/ electrons. What is physically relevant is an 
intermediate picture, where the fast non-/ electrons are 
allowed to screen the charge hole, while the /-electrons 
are not allowed to respond to the removed charge, a sit- 
uation which is not easily implemented in a total energy 
scheme. 

According to Janak's theorem^ the derivative of the 
LSD energy functional with respect to orbital occupation 
is equal to the LSD eigenvalue for the corresponding or- 
bital, which can be taken as a good representative of the 
removal energy of an electron in that state, if the charge 
of the orbital is everywhere small, as it will be for a band- 
like state. A similar theorem holds in SIC-LSD^ where 
the SIC eigenvalue 

e S ic=<f\H LSD + V SI c\f> (13) 

represents the derivative of the SIC-LSD total energy 
with respect to the occupation of the 'canonical orbitals', 
which in the context of a periodic solid means a Bloch 
state formed from the localized /-states on all the rare- 
earth sites. Here, Vsic represents the aditional poten- 
tial term arising in the self-consistency equation due to 
the self- interaction term in Eq. (@J. Hence, using the 
SIC eigenvalues, the occupied /-states (without multi- 
plets) appear as sharp resonances (SIC-LSD 'bands'—) 
below the valence bands (at around —11.5 eV for the 
case of Eu in GaN depicted in Fig. |3J, which is unrealis- 
tically deep for comparison to physical removal energies 
as observed in photoemission. The reason is that pho- 
toemission knocks out a localized /-electron rather than 
a 'canonical' /-state. 

Similar to the previous subsection, we may however 
apply a physically more reasonable but theoretically less 
rigorous approach by placing the removal energies of /- 
states at the transition state position. The eigenvalue 
given by Eq. (fTTTft represents the energy cost due to 
removal of the first infinitesimal part of the /-electron 
(since Vsic is evaluated for the initial ground state, i.e. 
at full occupancy), but as more and more of the /- 
electron is removed, the SIC potential term decreases, 
and eventually, only Hlsd is left. As an average we may 
therefore take the / removal energy as midway between 
the SIC-LSD and LSD band positions: 

e TS = i(< f\H LSD + Vswlf > + < f\H LSD \f >)• 

(14) 

In effect, the SIC potential is only counted with half its 
strength in the transition state approximation to the re- 
moval energy. By evaluating Hlsd in the initial state, 
i.e. without the hole in the / shell, we avoid the afore- 
mentioned effect of the /-hole pulling the / levels down. 
The transition state philosophy was also implemented in 




Ce Pr Nd Pm Sm Eu Gd Tb Dy Ho Er Tm Yb 



FIG. 10: (Color online) Transition state evaluation of the /- 
removal energies in GaAs (diamonds) and GaN (solid circles). 
The units are eV and relative to the VBM. 



Ref. l5ll albeit in a different manner, by invoking the av- 
eraging factor of h already in the total energy functional, 
while we do it here only for the removal energy, Eq. I|14fl , 
after self-consistency. 

Figure ITU1 displays the trends of e^ s through the lan- 
thanide series in both GaAs and GaN. The gradual in- 
creasing binding of the /-shell is clearly reflected in the 
transition state, only broken at Gd in the middle of the 
series. Furthermore, the / removal energy is consistently 
smaller (less negative) in GaN compared to GaAs. A few 
experimental data corroborate the interpretation of the 
transition state as a good measure of the / removal en- 
ergy: In ref. l33l the photoemission spectra of Eu doped 
GaN reveal the signature of the / 6 — > / 5 transition in 
the interval between 5 and 10 eV below the Fermi edge, 
i.e. the lowest f 5 ( 6 H) multiplet is situated at -5 eV, in 
excellent agreement with the value -5.1 eV in Fig. ^| 
For GdN films, the Gd /-emission peaks at -7.8 eV 52 , 
compared to the transition state energy of dilute Gd in 
GaN of -6.2 eV in Fig. EU For YbN, the similar onset 
of /-emission is estimated to occur around -6 eVj^ com- 
pared to the value -7.1 in Fig. ^| In CeAs, /-emission 
occurs around -2.5 eV^ as compared to the transition 
state estimate of -2.1 eV for dilute Ce in GaAs. In GdAs, 
/-emission occurs at -9 eV^ while the estimate of Fig. 
El for Gd in GaAs is -8.9 eV. For ErAs, the /-emission 
starts at -5 eV^ while the transition state position of 
dilute Er in GaAs is calculated as -6.7 eV. 



IV. SUMMARY 

Trends in the electronic structure of rare-earth substi- 
tutional impurities in GaAs and cubic GaN have been 
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discussed based on self-interaction corrected local spin 
density calculations. The trivalent configuration is found 
as the ground state of most dopants, while the divalent 
acceptor level is found above the theoretical band gap, 
and within the experimental gap only for Eu, Tm and 
Yb in GaAs and for Eu in GaN. The magnetic interaction 
of the rare earth ion with the host states at the valence 
and conduction band edges has been investigated and 
found to be relatively weak in comparison with 3d impu- 
rities. Hence it is unlikely that the rare-earth dopant by 
itself may induce room-temperature ferromagnetism and 
gigantic magnetic enhancement, as observed in certain 
experiments.- 1 ' 4 To explain these experiments it seems 
imperative to include the interaction of the rare-earth 



dopants with other defects, native or external, which will 
be pursued in our future research. 
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